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Abstract 

Based on the semi-classical extended Thomas-Fermi approach, we study the mass dependence 
of the symmetry energy coefficients of finite nuclei for 36 different Skyrme forces. The reference 
densities of both light and heavy nuclei are obtained. Eight models based on nuclear liquid drop 
concept and the Skyrme force SkM* suggest the symmetry energy coefficient a sym = 22.90 ± 0.15 
MeV at A = 260, and the corresponding reference density is pa — 0.1 fm~ 3 at this mass region. The 
standard Skyrme energy density functionals give negative values for the coefficient of the I 4 term 
in the binding energy formula, whereas the latest Weizsacker-Skyrme formula and the experimental 
data suggest positive values for the coefficient. 
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I. INTRODUCTION 


Nuclear symmetry energy has attracted a lot attention in recent years. In addition to its 
importance in the study of nuclear physics such as nuclear exotic structures and reactions 
induced by unstable nuclei, the symmetry energy also plays a role in the study of nuclear 
astrophysics, such as the r-process and the properties of neutron stars [1]. Although a 
great effort has been devoted in recent decades to investigate the symmetry energy [2-8], 
the density dependence of the symmetry energy, even at sub-saturation density region, is 
not very well constrained. One usually obtains the information of the symmetry energy 
from nuclear dynamical behavior in reactions [9-11] and the static properties of finite nuclei 
such as nuclear masses [12-19] and neutron skin thickness [20-23], or from the asymmetric 
nuclear matter based on various effective interactions [24-26], or from the observables in 
nuclear astrophysics. To understand the behavior of symmetry energy in nuclear masses and 
nuclear matter, it is crucial to establish a reliable connection between the mass dependence 
of the symmetry energy coefficient a sym (A) of finite nuclei and the density dependence of 
the symmetry energy E sym (p ) of nuclear matter, with which the obtained symmetry energy 
coefficients from various macroscopic-microscopic or liquid-drop models could be used to 
constrain the behavior of E syja (p) at sub-saturation densities. 

It is known that the density functional theory is widely used in the study of the nuclear 
ground state which provides us with a useful balance between accuracy and computation 
cost, allowing large systems with a simple self-consistent manner. With the same energy 
density functional, both the E sym (p ) and a sym (A) could be self-consistently investigated. 
In the framework of the effective Skyrme interaction [27] and the extended Thomas-Fermi 
(ETF) approximation, the Skyrme energy density functional approach was proposed for the 
study of nuclear structure and reactions [28-30]. One of the advantages in this approach 
is that the Coulomb energy, the Wigner energy and the microscopic shell corrections can 
be ’’cleanly” removed from the binding energies of nuclei, which is important to accurately 
obtain the information of symmetry energy. Another advantage of this approach is that the 
higher-order terms of the symmetry energy can be investigated simultaneously. 

On the other hand, the investigation of symmetry energy is helpful to improve and test 
the reliability of various nuclear mass models, especially for the predictions of the masses 
of nuclei approaching neutron drip line. Both the macroscopic-microscopic mass formulas 
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[12, 13] and the microscopic Hartree-Fock-Bogoliubov (HFB) models [14, 15] can reproduce 
the measured masses of more than 2000 nuclei with an rms error of several hundred keV. 
However, the predictions for extremely neutron-rich nuclei are quite different from these 
models although the similar symmetry energy coefficient is adopted. For example, adopting 
almost the same symmetry energy coefficient J ~ 30 MeV, the predicted mass of 168 Sn with 
the WS4 model [13] is larger than those with the HFB17 model [14] by 24.4 MeV, which is 
difficult to be explained from the microscopic shell and pairing corrections. The higher-order 
term of symmetry energy might play a relevant role to the large difference of the predicted 
masses, in addition to other effects. It is therefore necessary to investigate the symmetry 
energy coefficient of finite nuclei especially the coefficient of the J 4 term in the standard 
Skyrme force. 


II. THE THEORETICAL MODEL 


The nuclear energy part of a nucleus at its ground state can be expressed as the integral 
of the Skyrme energy density functional 'H(r): 


E= ft[p q (r),T q (r), J q (r)] dr, 


( 1 ) 


with the local nucleon densities p q (r), kinetic energy densities r q (r) and spin-orbit, densities 
Jq( r ) (q = n for neutrons and q = p for protons). By using the extended Thomas-Fermi 
approach up to the second order in h (ETF2) [28-30], the kinetic energy density and the 
spin-orbit density can be expressed as a functional of nuclear density and its gradients. In 
the semi-classical ETF2 approach, one can obtain the ’’macroscopic” part of the nuclear 
energy E\p n (r), p p (r)\ of a nucleus. The microscopic shell, pairing, and Wigner effects are 
not involved in this semi-classical approach, which is helpful to study the symmetry energy 
coefficients of finite nuclei as cleanly as possible. 

To remove the influence of nuclear deformations, we use the spherical Fermi function 


Pg( r ) = 


$ 


( 2 ) 


l + exp(^) 

for describing the density of a nucleus. p^\ R q and a q denote the central density, radius, 
and surface diffuseness of nuclei, respectively. The central density is determined from the 
conservation of particle number. By using the optimization algorithm [30] and varying the 
four variables R p , a p , R n , a n in Eq.(2) for a given nucleus, one can self-consistently obtain 
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FIG. 1. (Color online) Binding energy per particle as a function of I = (N — Z)/A. The Coulomb 
interaction is not involved in all calculations with the ETF2 approach. The squares denote the 
results with the Skyrrne force SkM*, and the curves denote the fit with Eq.(3). 

the minimal ’’macroscopic” energy E of the nucleus. It is found that the ETF approach 
can reproduce quite accurately the Hartree-Fock average field and corresponding density 
distribution of 208 Pb [28]. As the same as those did in Ref. [31], the calculations have been 
carried out for nuclei with huge numbers of nucleons, of the order of 10 6 , in order to perform 
a reliable extrapolation in the inverse radius, and the Coulomb interaction has been ignored 
to be able to approach nuclei of arbitrary sizes and to avoid radial instabilities characteristic 
of systems with very large atomic numbers. 

III. RESULTS AND DISCUSSIONS 

As an example, ’’macroscopic” energy per particle of nuclei is studied by adopting the 
Skyrrne force SkM* [32] together with the ETF2 approach, since SkM* is very successful for 
describing the bulk properties and surface properties of nuclei. Figure 1 shows the calculated 
’’macroscopic” energy per particle E/A as a function of isospin asymmetry I = (N — Z)/A 
for nuclei with given mass number A. The calculations are performed in the region of 
—0.2 < I < 0.3 which generally covers all nuclei with known masses. The squares denote 
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FIG. 2. (Color online) Symmetry energy coefficients as a function of mass number of nuclei. 

the calculated results with the Skyrme force SkM*. Ignoring the Coulomb energy and 
microscopic corrections, the binding energy per particle of nuclei is usually written as 

E/A = eo(A) + a sym (A)I 2 + ai%(A)I 4 + (?(/ 6 ). (3) 

The solid curves in Fig. 1 denote the fit to the squares by using the form of Eq.(3). The 
results of SkM* can be remarkably well reproduced by the solid curves, with the rms devi¬ 
ations smaller than 10” 6 MeV. Figure 2 shows the obtained symmetry energy coefficients 
of finite nuclei as a function of mass number. The squares and circles denote the extracted 
results from the Skyrme force SkM* for the values of a sym and ai yL, respectively. We find 
that the values of oiyL are negative from the calculations of SkM* and the values of a sym 
approach 30 MeV with increasing of mass number. To check the convergence of Eq.(3), we 
re-calculate the values of a sym by neglecting the J 4 and higher-order terms in Eq.(3) and 
End that the obtained values of a sym are slightly reduced by about 0.2 MeV. 

To further check the asymptotic behavior of the symmetry energy coefficient when A —y 
oo, we further use a Polynomial fit to reproduce the calculated a sym and ai y m . The solid 
curves in Fig. 2 denote the results of the Polynomial fit to the scattered symbols. To see the 
results more clearly, the results for a sym with SkM* (open squares) and with SLy4 [33] (open 
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FIG. 3. (Color online) Symmetry energy coefficient a sym as a function of A _1// ' 3 . The open circles 
and squares denote the results of ETF2 plus SLy4 and SkM*, respectively. The curves denote the 
Fit with Eq.(4). 

circles) are shown in Fig. 3 as a function of A~ 1 / 3 . According to the liquid drop model, the 
symmetry energy coefficient of a finite nucleus is usually written as 

®sym(^4) = J — o ss A 1 ^ 3 + a cs A 2 / 3 (4) 

by using the Leptodermous expansion in terms of powers of A^ 1 ^ 3 . J ^ 28 — 34 MeV denotes 
the symmetry energy of nuclear matter at normal density. a ss is the coefficient of the surface- 
symmetry term and a cs denotes the curvature-symmetry term. Through the Polynomial fit 
to the results of different Skyrme forces, one can obtain the corresponding values of J, a ss 
and a cs . 

For nuclear matter, the symmetry energy in the standard Skyrme energy density func¬ 
tional is expressed as 

, . 1 h 2 f 3 n 2 \ 2 ^ 3 2/3 1 . . 

E ‘* M = 3 2^ {—) P 8‘ o(2l ° + 1)P 

1 / 37T 2 \ 2 ^ 3 1 

- 24 \ ~ 2 J 0symp5/3 ~~ 48 t3 ^ 2x3 + ^ pa+1 
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with 0 sym = 3tiXi —^(4 + 5 x 2 ). to, tii ^ 2 , t 3 , xo, x\, X 2 , X 3 and a are the Skyrme parameters. 
For the parameter set SkM* and SLy4, the corresponding symmetry energy E syjn (p 0 ) at 
saturation density p 0 from Eq.(5) is 30.0 and 32.0 MeV, respectively. From Fig. 3, we also 
note that the obtained value of J from the symmetry energy coefficient of finite nuclei in 
the ETF2 approach is very close to the corresponding value of E sym (p 0 ). Furthermore, the 
obtained binding energies per particle eo(oo) from the extrapolation with the Polynomial 
fit, i.e., the lowest energy in Fig. 1 when the mass number A — y 00 , are —15.78 MeV for 
SkM* and —15.99 MeV for SLy4, which are also very close to the corresponding energy per 
nucleon E/A(p 0 ) at saturation density of asymmetric nuclear matter [34]. 

In this work, we have performed a systematical study of the symmetry energy coefficient 
for 36 different Skyrme forces in which the corresponding incompressibility coefficient for 
symmetry nuclear matter is about = 240 ± 30 MeV and the measured masses of 54 
spherical nuclei according to the predicted shapes of nuclei from the WS4 model can be 
roughly reproduced by using the ETF2 approach (with the rms deviations to the masses 
of the 54 nuclei less than 10 MeV). The results are listed in Table I. We find that the 
average deviation between J obtained from the symmetry energy coefficient of finite nuclei 
and E sym (p 0 ) from Eq.(5) is only 0.04 MeV. In addition, it is important and necessary to 
check the obtained coefficient ai y L of the J 4 term from the symmetry energy coefficients of 
finite nuclei. The fourth-order term of the energy per nucleon at saturation density E satj 4 
of asymmetric nuclear matter were systematically investigated by Lie-Wen Chen et al. for 
different Skyrme forces [48]. With the Polynomial fit to the calculated ai y L for a certain 
Skyrme force, one can obtain the corresponding value of oiyL(oo) which is also listed in 
Table I. We find that aiy m (oo) obtained in the ETF2 approach is close to the corresponding 
value of E sat p. The average deviation between a| y m(oo) extracted in this work and C sat; 4 
given in Ref. [48] is only 0.23 MeV. These tests for J, eo(oo) and aiy m (oo) indicate that the 
ETF2 approach is reliable for the study of the symmetry energy coefficient of finite nuclei. 

With the ETF2 approach, the connection between the symmetry energy coefficient 
a sym (7l) and the symmetry energy of nuclear matter E sym (p) can be established by using the 
relationship a sym (A) = E sym (pA) [ 20 ], where pa is the reference density. It is usually thought 
that the reference density of a heavy nucleus such as 208 Pb is about 0.1 firm 3 [17, 20]. In Fig. 
4, we show the calculated reference density pa with the ETF2 approach by adopting the pa¬ 
rameter set SkM*. According to the results of SkM*, we note that the reference density for 
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TABLE I. Properties of nuclear symmetry energy with the ETF2 approach (in MeV). 


Force 

L c 

J 

^SS 

^CS 

&sym(oo) 

agL(208) 

aiyL(40) 

BSkl [35] 

30.31 

27.95 

28.9 

7.3 

0.3 

-0.8 

-1.4 

MSk2 [36] 

39.66 

29.99 

40.0 

19.2 

-2.2 

-2.6 

-2.6 

MSk6 [36] 

30.92 

28.01 

29.3 

7.8 

0.2 

-0.8 

-1.4 

RATP [37] 

38.06 

29.25 

40.5 

20.2 

-1.4 

-2.3 

-2.4 

SkM [38] 

44.55 

30.71 

50.1 

32.9 

-5.2 

-3.5 

-2.6 

SkM* [32] 

42.93 

30.00 

50.2 

33.0 

-4.4 

-2.9 

-2.2 

SkMP [39] 

49.64 

29.79 

53.3 

40.2 

-10.0 

-5.2 

-3.3 

SkSCl [40] 

28.87 

28.11 

26.7 

4.6 

0.5 

-0.4 

-1.2 

SkSC4 [41] 

29.06 

28.81 

27.0 

4.0 

0.5 

-0.5 

-1.4 

SkTl* [42] 

47.69 

31.97 

56.0 

39.4 

-6.6 

-4.7 

-3.4 

SkT3 [42] 

46.80 

31.46 

48.7 

29.7 

-6.4 

-5.4 

-4.4 

SkT3* [42] 

47.25 

31.64 

52.6 

35.0 

-6.6 

-5.4 

-4.3 

SkT4 [42] 

61.79 

35.24 

78.1 

71.0 

-18.3 

-8.5 

-4.8 

SkT6 [42] 

38.85 

29.96 

41.2 

20.1 

-1.8 

-2.4 

-2.5 

SkT7 [42] 

38.68 

29.51 

40.1 

19.7 

-1.5 

-2.4 

-2.6 

SkT8 [42] 

38.81 

29.92 

41.4 

20.8 

-2.1 

-2.5 

-2.5 

Skz2 [43] 

35.24 

32.02 

37.5 

13.6 

0.2 

-1.2 

-1.9 

Skz3 [43] 

32.22 

32.03 

31.8 

6.6 

0.1 

-1.6 

-2.6 

Skz4 [43] 

28.51 

32.02 

26.0 

0.3 

0.1 

-1.7 

-3.0 

SLyO [44] 

41.90 

31.97 

49.3 

29.1 

-4.7 

-3.5 

-2.8 

SLyl [44] 

41.89 

31.97 

51.4 

32.0 

-4.7 

-3.0 

-2.2 

SLy2 [44] 

41.99 

31.99 

47.4 

26.5 

-4.7 

-4.1 

-3.5 

SLy3 [44] 

41.34 

31.98 

50.0 

30.1 

-4.3 

-3.0 

-2.3 

SLy4 [33] 

41.59 

31.99 

50.4 

30.7 

-4.5 

-3.0 

-2.3 

SLy5 [33] 

42.19 

31.98 

51.6 

32.2 

-4.9 

-3.1 

-2.2 

SLy6 [33] 

41.96 

31.94 

47.6 

27.6 

-4.8 

-3.4 

-2.7 

SLy7 [33] 

41.86 

31.98 

46.6 

26.2 

-4.7 

-3.5 

-2.8 






TABLE II. ( Continued .) 


Force 

L c 

J 

^SS 

^CS 

&sym(oo) 

Osym(208) 

Osym (40) 

SLy8 [44] 

41.89 

31.98 

50.6 

31.1 

-4.7 

-3.1 

-2.3 

SLy9 [44] 

44.95 

31.95 

51.0 

32.7 

-6.4 

-4.3 

-3.2 

SLylO [33] 

39.39 

31.98 

41.4 

19.7 

-3.1 

-3.4 

-3.3 

SLy230a [45] 

39.39 

31.96 

45.2 

23.4 

-4.4 

-3.5 

-3.0 

SLy230b [45] 

41.60 

31.99 

50.4 

30.7 

-4.4 

-3.0 

-2.3 

SV-sym32 [46] 

49.29 

32.12 

50.5 

32.5 

-6.7 

-6.0 

-5.0 

V080 [47] 

30.12 

28.01 

29.8 

8.7 

0.8 

-0.2 

-0.9 

V090 [47] 

30.31 

28.01 

29.0 

7.6 

0.6 

-0.5 

-1.2 

V110 [47] 

30.15 

28.01 

28.7 

7.1 

0.3 

-0.7 

-1.4 



40 80 120 160 200 240 280 


Mass number A 


FIG. 4. Reference density obtained from SkM* as a function of mass number of nuclei, 
nuclei with A = 260 is about 0.1 fm ~ 3 . In Table I, we also list the slope L c = 3 p c ( dE * yrn "j 

V P J p=pc 

of the symmetry energy E sym at the sub-saturation density of p c = 0.1 firm 3 . We note that 
the value of the surface-symmetry coefficient a ss increases linearly with the value of L c in 
general, with which the corresponding value of L c for different macroscopic-microscopic and 
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FIG. 5. (Color online) Symmetry energy coefficients of nuclei extracted from nuclear masses, 
liquid-drop models could be estimated. 

To illustrate the importance of the slope L c and the density p c , we simultaneously inves¬ 
tigate the symmetry energy coefficients a sym (H) extracted by using various liquid drop and 
macroscopic-microscopic mass models together with the measured masses of nuclei in this 
work. It is known that only one-third of the nucleons in heavy nuclei occupy the saturation 
density area [49]. Consequently, nuclear observables related to the average properties of 
nuclei, such as masses or radii, constrain the equation of state not at the saturation density 
but rather around the so-called ’’critial” density p c ~ 0.1 fm -3 . In Refs. [17-19], the authors 
determine the coefficients of the symmetry energy term in the liquid drop model, assuming 
the form a sym (A) = J — a ss H -1 / 3 or a sym (H) = J/(l + k/R 1 / 3 ) [50] and subtracting the 
Coulomb and Wigner terms from the measured binding energies of nuclei. The obtained 
symmetry energy coefficients a aym (H) are shown in Fig. 5. Here, we also show the corre¬ 
sponding a sym adopted in five different macroscopic-microscopic mass models [13, 51? -53] 
for comparison. Because the nuclear symmetry energy in nuclear masses is highly correlated 
with the Coulomb energy, the Wigner energy, and as well as the shell corrections, one can 
see from Fig. 5 that the extracted symmetry energy coefficients are quite different for light 
nuclei and the nuclei with huge numbers of nucleons. However, for nuclei with A = 260, 
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FIG. 6. (Color online) Binding energy per particle after removing the Coulomb energy, Wigner 
energy and traditional symmetry energy. 

all these different models give quite similar predictions a sym = 22.90 ± 0.15 MeV. With the 
ETF2 approach, the corresponding result from SkM* is a sym = 22.95 MeV, which is in good 
agreement with the results in Fig. 5. 

As mentioned previously, the predicted masses of extremely neutron-rich nuclei are quite 
different from the WS4 and HFB17 models. To understand the large deviations, we analyze 
the influence of the / 4 term on the masses of neutron-rich nuclei. According to the well 
known liquid-drop formula, the ground state energy of a nucleus is expressed as 

Eb ~ eo A + Ec + E w + a sym I 2 A + a^^I^A + ..., (6) 

with eo = a„ + a s A -1 / 3 for the volume and surface terms, the Coulomb energy Ec = 
0.71Z 2 /A 1//3 (1 — 0.76Z -2 / 3 ) and the Wigner energy Ew = 47|/| [54], After removing the 
Coulomb energy, the Wigner energy, and the traditional symmetry energy ( a sym I 2 A ) from 
the total energy Eb of a nucleus, one obtains the binding energy per particle of finite nuclei 

£ 0 ~ e 0 + dgyln/ 4 , (7) 

neglecting the microscopic corrections in nuclei. In Fig. 6 , we show the calculated £q for 
nuclei with A = 160. The choice of A = 160 is to avoid the influence of the shell effect as 
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possible. The black squares denote the experimental data. The circles, short-dashed and 
solid curves denote the predictions of WS4 [13], WS* [12] and HFB17 [14], respectively. 
We take a sym (160) = 21.7 MeV in the calculations according to the predictions of the eight 
models in Fig. 5. In addition, we show the results (dashed curve) of ETF2 together with 
SkM* for comparison. Here, the results of SkM* is shifted by 0.17 MeV to reproduce the 
experimental data. One can see that for nuclei with known masses, all the three mass models 
WS4, WS* and HFB17 can reproduce the experimental data remarkably well. However, the 
trend becomes quite different for extremely neutron-rich nuclei. The calculated e 0 from the 
HFB17 model decreases with the isospin asymmetry, whereas the result of WS4 increases 
with the isospin asymmetry. The trend of £$ from the Skyrme force SkM* is similar to 
those from HFB17 which is based on the Skyrme force Bskl7 [14]. With the help of the 
ETF2 approach, we obtain the coefficient of the J 4 term ai y L = —2.84 MeV for nuclei with 
A = 160 by adopting SkM*. The corresponding values of the coefficient ai y ]n for nuclei 
with mass number A = 208 and A — 40 are also listed in Table I, denoted by ag } m (208) 
and ai y m(40), respectively. For almost all 36 selected Skyrme forces, the coefficients ai y L for 
finite nuclei are negative. In the WS4 model, the surface diffuseness correction is taken into 
account for unstable nuclei, which causes the enhancement of the symmetry energy for nuclei 
approaching drip lines and thus leads to the enhancement of £q at large isospin asymmetry 
region. In the WS* model, neither the surface diffuseness correction nor the I 4 term is 
considered. Therefore, the results of e 0 from WS* look flat in general. The decrease of the 
value of £ 0 at small isospin asymmetry (/ < 0.1) in the macroscopic-microscopic and HFB17 
calculations could be due to the influence of shell effects around magic number N = 82. 

Very recently, Jiang et al. [55] extracted the coefficient of the J 4 term with the double 
difference of the symmetry energy term together with the measured masses (AME2012) [56], 
and a value of ai y L = 3.28 MeV was obtained. Here, we also show in Fig. 6 the trend of 
e 0 with a s y L = 3.28 MeV (dot-dashed curve). One sees that the trend of £ 0 from Jiang et 
al. are close to those from the WS4 model with which the rms deviation to 2353 measured 
masses [56] is only 298 keV. With the same approach proposed in Ref. [55], we obtain 
Osym = —3.07 MeV for the HFB17 mass model, which is generally consistent with the result 
of the ETF2 approach. For 168 Sn (I ~ 0.4) mentioned previously, the contribution of the 
I 4 term to the binding energy is aiylnI 4 A ps ±13 MeV according to the values of ai ym in 
HFB17 and WS4. It seems that the opposite values for the coefficient of / 4 term used in 
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the HFB17 and WS4 models result in the large difference at the predictions of the masses 
of heavy nuclei approaching the neutron drip-line. 

IV. SUMMARY 

Based on the extended Thomas-Fermi (ETF2) approximation together with the restricted 
density variational method, the symmetry energy coefficients of finite nuclei including the 
coefficients of the I 4 term have been systematically investigated with 36 different Skyrme 
energy density functionals. From nuclei with mass number A = 20 to the nuclei with huge 
numbers of nucleons, of the order of 10 6 , we study the mass dependence of the symmetry 
energy coefficients. With the extrapolations from the calculated results for finite nuclei, the 
asymptotic values of the symmetry energy coefficient J, of the binding energy per particle 
eo(oo) and of the coefficient OsyL(oo) in the I 4 term when A —» oo have been compared with 
the corresponding values from the asymmetric nuclear matter. The obtained results from 
the finite nuclei are very close to those from the nuclear matter, which indicates that the 
ETF2 approach is reliable for extracting the symmetry energy coefficients of finite nuclei, 
because the Coulomb energies, the Wigner energies and the shell corrections can be ’’cleanly” 
removed in the calculations. 

With the help of the ETF2 approach, the reference density of finite nuclei is also in¬ 
vestigated. We find that the extracted symmetry energy coefficients a sym (A) from different 
liquid-drop and macroscopic-microscopic models for heavy nuclei such as nuclei with A = 260 
are in good agreement with each other, a sym (260) = 22.90 ± 0.15 MeV, although the un¬ 
certainty for very light nuclei and those for nuclear matter is relatively large. The ETF2 
approach with SkM* gives a similar result a sym (260) = 22.95 MeV and the corresponding 
reference density is pa — 0.1 fm -3 . 

In addition, we analyze the large deviation between the HFB17 and WS4 model for the 
predictions of the masses of extremely neutron-rich nuclei. From the predicted binding 
energy per particle for nuclei with A = 160 after removing the Coulomb energy, Wigner 
energy and traditional symmetry energy, we find that the HFB17 model gives a negative 
value for the coefficient ai y L of the / 4 term, whereas the WS4 model gives a positive value 
due to considering the surface diffuseness correction of nuclei. The extracted result from 
nuclear masses with the double difference to the symmetry energy term suggests a positive 
value for ai ym . The negative values for ai y L from the systematic study of all 36 different 
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Skyrme forces imply that the experimental data used in the different fitting protocols do not 
constrain this part of the standard Skyrme functionals and/or new terms should be included 
in future Skyrme functionals. 


ACKNOWLEDGEMENTS 

This work was supported by National Natural Science Foundation of China (Nos 
11275052, 11365005, 11422548, 11475004, 11225524, 11305101) and the 973 Program of 
China (Grant Nos. 2013CB834401). N. W. is grateful to Lie-Wen Chen and Li Ou for 
valuable discussions, and acknowledges the support of the Open Project Program of State 
Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy 
of Sciences, China (No. Y4KF041CJ1). 


[1] B. A. Li, A. Ramos, G. Verde, I. Vidaha, ” Topical issue on nuclear symmetry energy ”, Eur. 
Phys. J. A50, 9 (2014); and references therein. 

[2] B. A. Li, L.W. Chen, and C. M. Ko, Phys. Rep. 464, 113 (2008). 

[3] L. W. Chen, C. M. Ko, and B. A. Li, Phys. Rev. Lett. 94, 032701 (2005). 

[4] D. V. Shetty, S. J. Yennello, and G. A. Souliotis, Phys. Rev. C 76,024606 (2007). 

[5] A. S. Botvina, O. V. Lozhkin, and W. Trautmann, Phys. Rev. C 65,044610 (2002). 

[6] A. W. Steiner and S. Gandolfi, Phys. Rev. Lett. 108 , 081102 (2012). 

[7] A. W. Steiner, M. Prakash, J. M. Lattimer, P. J. Ellis, Phys. Rep. 411, 325 (2005). 

[8] J. Dong, W. Zuo, and W. Scheid, Phys. Rev. Lett. 107 , 012501 (2011). 

[9] Y. X. Zhang, P. Danielewicz, M. Famiano, Z. Li, W. G. Lynch, and M. B. Tsang, Phys. Lett. 
B 664, 145 (2008). 

[10] M. B. Tsang, Y. X. Zhang, P. Danielewicz, M. Famiano, Z. Li, W. G. Lynch, and A. W. 
Steiner, Phys. Rev. Lett. 102 , 122701 (2009). 

[11] L. Trippa, Gianluca Colo, and E. Vigezzi, Phys. Rev. C 77, 061304(R) (2008). 

[12] N. Wang, Z. Liang, M. Liu and X. Wu, Phys. Rev. C 82 , 044304 (2010). 

[13] N. Wang, M. Liu, X. Z. Wu, and J. Meng, Phys. Lett. B 734, 215 (2014). 

[14] S. Goriely, N. Charnel and J. M. Pearson, Phys. Rev. Lett. 102 , 152503 (2009). 

[15] S. Goriely, N. Charnel, and J. M. Pearson, Phys. Rev. C 88, 061302(R) (2013) 


14 



[16] P. W. Zhao, Z. P. Li, J. M. Yao, and J. Meng, Phys. Rev. C 82, 054319 (2010). 

[17] M. Liu, N. Wang, Z.-X. Li, and F.-S. Zhang, Phys. Rev. C 82 , 064306(2010). 

[18] H. Jiang, G. J. Fu, Y. M. Zhao, and A. Arima, Phys. Rev. C 85 , 024301 (2012). 

[19] J. L. Tian, H. T. Cui, K. K. Zheng, and N. Wang, Phys. Rev. C 90, 024313 (2014). 

[20] M. Centelles, X. Roca-Maza, X. Vinas, and M. Warda, Phys. Rev. Lett. 102 , 122502 (2009). 

[21] N. Wang, L. Ou, and M. Liu, Phys. Rev. C 87, 034327 (2013). 

[22] N. Wang and T. Li, Phys. Rev. C 88, 011301(R) (2013). 

[23] Z. Zhang, L. W. Chen, Phys. Lett. B 726, 234 (2013). 

[24] W. Zuol, I. Bombaci, and U. Lombardo, Eur. Phys. J. A 50 , 12 (2014). 

[25] D. Vretenar, T. Niksic, and P. Ring, Phys. Rev. C 68, 024310 (2003). 

[26] V. Greco, V. Baran, M. Colonna, M. Di Toro, T. Gaitanos, H.H. Wolter, Phys. Lett. B 562, 
215 (2003) 

[27] D. Vautherin, D. M. Brink, Phys. Rev. C 5, 626 (1972). 

[28] J. Bartel, K. Bencheikh. Eur. Phys. J 14 , 179 (2002). 

[29] M. Brack, C. Guet, and H.-B. Hakanson, Phys. Rep. 123 , 275 (1985). 

[30] M. Liu, N. Wang, Z. Li, X. Wu, and E. Zhao, Nucl. Phys. A 768, 80 (2006). 

[31] P.-G. Reinhard, M. Bender, W. Nazarewicz, and T. Vertse, Phys. Rev. C 73 , 014309 (2006). 

[32] J. Bartel, P. Quentin,M. Brack, C. Guet, and H.-B. Hakansson, Nucl. Phys. A 386 , 79 (1982). 

[33] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A 635 , 231 
(1998). 

[34] M. Dutra, O. Lourengo, J. S. Sa Martins, and A. Delfino, et al., Phys. Rev. C 85 , 035201 

( 2012 ). 

[35] M. Samyn et al., Nucl. Phys. A 700 , 142 (2002). 

[36] F. Tondeur, S. Goriely, J. M. Pearson, and M. Onsi, Phys. Rev. C 62 , 024308 (2000). 

[37] M. Rayet, M. Arnould, F. Tondeur, and G. Paulus, Astron. Astrophys. 116 , 183 (1982). 

[38] H. Krivine, J. Treiner, and O. Bohigas, Nucl. Phys. A 336 , 155 (1980). 

[39] L. Bennour, P.-H. Heenen, P. Bonche, J. Dobaczewski, and H. Flocard, Phys. Rev. C 40 , 2834 
(1989). 

[40] J. M. Pearson, Y. Aboussir, A. K. Dutta, R. C. Nayak, M. Farine, and F. Tondeur, Nucl. 
Phys. A 528, 1 (1991). 

[41] Y. Aboussir, J. M. Pearson, A. K. Dutta, and F. Tondeur, Nucl. Phys. A 549, 155 (1992). 


15 


[42] F. Tondeur, M. Brack, M. Farine, and J. M. Pearson, Nucl. Phys. A 420 , 297 (1984). 

[43] J. Margueron, J. Navarro, and N. Van Giai, Phys. Rev. C 66, 014303 (2002). 

[44] E. Chabanat, Ph.D. Thesis, University of Lyon, 1995. 

[45] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A 627 , 710 
(1997). 

[46] P. Kliipfel, P.-G. Reinhard, T. J. Btirvenich, and J. A. Maruhn, Phys. Rev. C 79 , 034310 
(2009). 

[47] J. M. Pearson and S. Goriely, Phys. Rev. C 64 , 027301 (2001). 

[48] L. W. Chen, B. J. Cai, C. M. Ko, B. A. Li, C. Shen, and J. Xu, Phys. Rev. C 80, 014322 
(2009). 

[49] E. Khan, J. Margueron, and I. Vidana, Phys. Rev. Lett. 109 , 092501 (2012). 

[50] P. Danielewicz and J. Lee, Nucl. Phys. A 818 , 36 (2009). 

[51] H. F. Zhang, J. M. Dong, N. N. Ma, et ah, Nucl. Phys. A 929 , 38 (2014) 

[52] M. W. Kirson, Nucl. Phys. A 798, 29 (2008). 

[53] A. Bhagwat, X. Vinas, M. Centelles, P. Schuck, and R. Wyss, Phys. Rev. C 81, 044321 (2010). 

[54] W. Satula, D. J. Dean, J. Gary, S. Mizutori, W. Nazarewicz, Phys. Lett. B 407 , 103 (1997). 

[55] H. Jiang, M. Bao, L. W. Chen, Y. M. Zhao, and A. Arima, Phys. Rev. C 90, 064303 (2014). 

[56] G. Audi, M. Wang, A. H. Wapstra, F. G. Kondev, M. Mac-Cormick, X. Xu, and B. Pfeiffer, 
Chin. Phys. C 36, 1287 (2012). 


16 


